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We present the exact analysis of a chiral Brownian motor and heat pump. Optimization of the 
construction predicts, for a nanoscale device, frequencies of the order of kHz and cooling rates of 
the order of femtojoule per second. 

PACS numbers: 05.70.Ln, 05.40.Jc, 07.10.Cm, 07.20. Pe 



Brownian motors have been studied intensively since 
the early 1990s [T]. This interest coincided with develop- 
ments in bioengineering and nanotechnology, where un- 
derstanding and designing a motor in the shape of a small 
biological or artificial device is an important issue. Most 
of the motors investigated in this context are powered by 
chemical energy. Brownian motors driven by a tempera- 
ture gradient [5] have a fundamental appeal, since their 
operation is directly related to basic questions such as 
Carnot efficiency. Maxwell demons and the foundations 
of statistical mechanics and thermodynamics [3|. The 
additional significance of the thermal Brownian motor 
comes from the recent observation that it can operate as 
a refrigerator |3], see also [5]- In fact, this property is, 
at least in the regime of linear response, a direct conse- 
quence of Onsager symmetry: if a temperature gradient 
generates motion, an applied force will generate a heat 
flux. This principle is well known in its application to 
electro-thermal devices, displaying the Peltier, Seebeck 
and Thompson effects At variance however with 

these macroscopic devices, rectiflcation of nonequilibrium 
thermal fluctuations provide the driving mechanism for 
Brownian refrigeration. The latter become more promi- 
nent, and so do the resulting motor and cooling functions, 
as the apparatus becomes smaller. 

Previous models for the Brownian refrigerator assume 
translational motion of the engine [J. This construction 
obviously poses difficulties in its technological implemen- 
tation, while the resulting friction is expected to lower 
the efficiency. In this letter, we present a chiral rota- 
tional model, in which these problems do not occur, and 
which has the extra benefit that the choice of the axis of 
rotation provides an additional parameter that can be op- 
timized. A related question, that will also be addressed, 
is the optimal chiral shape of the engine. The observed 
optimized rates of rotation and heat transfer are found 
to be significantly larger than in the translational coun- 
terpart, suggesting the technological implementation of 
such devices. 

Since the properties of the Brownian heat pump follow 
by Onsager symmetry from those of the Brownian mo- 
tor, we first focus on the latter. The basic construction is 
represented in Fig.[lja). The engine consists of two parts 
linked by a rigid axis (which we take to be the 2;-axis), 
around which the whole construction is free to rotate as 




FIG. 1: (Color online) (a) Schematic representation of the 
chiral motor (mass AI) rotating as a single unit along the ver- 
tical z-axis. Each motor element resides in a separate com- 
partment, filled with gas particles (mass m) at temperatures 
Ti, T2 and densities pi, p2, respectively, (b) The planar shape 
of one motor part is determined by the position vector f{x, y) 
of its perimeter (the origin being the axis of rotation). The 
collision rule Eq. ([2| is given in terms of the polar angle if, 
determining the direction orthogonal to the perimeter. The 
properties of the engine are expressed in terms of the tangen- 
tial component r„ — f-e„ of r along the perimeter cf. Eqs. ([s] 
[5]|8|. Maximal rotation and refrigeration speed is attained for 
one motor unit being the enantiomorph of the other. The cor- 
responding optimal shape (and axis implantation marked by a 
star) is the spiral form depicted here for maximum frequency 
(b) and maximum net cooling power (c). 



a single entity. This rotational motion is induced, follow- 
ing Newton's laws, by the collisions with surrounding gas 
particles. The question of interest is under which condi- 
tions sustained rotational motion will arise. Following 
the Curie principle, the breaking of symmetry plays a 
crucial role. There are two obvious symmetries involved: 
the statistical symmetry of the microscopic dynamics and 
the chiral symmetry of the device. If the motor units are 
achiral, in the sense that clockwise and counterclockwise 
rotation cannot be distinguished, no sustained motion 
will appear. When both motor units reside in a single 
compartment at equilibrium, sustained rotational motion 
will not appear, even for chiral motor units, because it 
would violate the second law of thermodynamics, or refer- 
ring to the basic underlying symmetry, because detailed 
balance should hold [7] . Both symmetries will be broken 
if we consider chiral units residing in separate compart- 
ments at unequal temperatures Ti . The index i runs over 
the different reservoirs, cf. Fig. [IT a) for a schematic rep- 
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resentation in the case of two reservoirs i = 1,2. As 
we proceed to show, the resulting average rotational fre- 
quency can be calculated exactly from microscopic dy- 
namics, at least in a limiting case. For simplicity, we will 
restrict the theoretical analysis in this letter to the case 
of a two-dimensional device. The corresponding results 
for three-dimensional cylindrical objects, as depicted in 
Fig.jlja), are obtained by appropriate rescaling with the 
hight h of the units (for more details, see [5]). 

Turning to the exact microscopic analysis, we con- 
sider convex (two-dimensional) units, residing in reser- 
voirs that are infinitely large and are filled with dilute 
gases at equilibrium. The reservoirs play the role of ideal 
thermostats with which the engine is exchanging energy. 
In the limit of high dilution and a heavy engine (mass 
M, moment of inertia /), the collisions of the engine 
parts with the gas particles (mass m) become uncorre- 
lated events and the following exact Boltzmann master 
equation for the probability distribution Pt{Lo) of its an- 
gular velocity holds: 



dt 



''du' [W^l^,Pt{L^') - W^,l^PtiL0)] . (1) 

Here Wf^\^^i is the transition probability per unit time for 
the motor to change its angular velocity from uo' to lo 
by one collision. We will assume that the collisions are 
perfectly elastic and that the interaction force is short- 
ranged and central. Conservation of the total energy and 
of the total angular momentum in the 0-dircction of the 
colliding pair, and of the tangential component of the 
momentum of the gas particle, leads to the following col- 
lision law: 

^ I , . (2) 

X sm ip — y cos ip + [x sm ip — y cos (p) 

specifying uj in terms of its pre-coUisional value lo' and 
the pre-coUisional speeds v' = {v'^, v'y) of the gas particle. 
(p is the polar angle of the surface at impact, see Fig.jljb). 

The transition probability W^^^i^r can now be calculated 
following standard arguments from the kinetic theory of 
gases. Taking into account that the velocity distributions 
of the particles 0i(v) are Maxwellian at the density pi 
and temperature Ti of their bath, one finds the following 
explicit expression: 



X exp 



({r, + —)ALu + 2r„Lu'^ / Uvl) 



(3) 



The contour integral f dlt is over the perimeter of the 
engine part in each reservoir, which can be of any convex 
shape. We have introduced the tangential component of 
the position vector r along the perimeter, r,, = f ■ e„ = 



—X sin tp + y cos p cf. Fig. [l|b), the thermal speed of the 
gas particles vt^ — ^/SkTi/nm, and the angular velocity 
increment Alu — lj — lj' . 

In equilibrium, Ti = T, one easily veri- 
fies that the Boltzmann distribution P'^'^{uj) = 
y^I/{2nkT) exp[-/w7(2/cr)] is the unique steady 
state solution of Eq. Q. As required by statistical 
mechanics, this solution satisfies detailed balance: 
W^l^,P''^iLo') = W_^,\_^P-^{~uo). 

A general solution away from equilibrium is not avail- 
able, hence we resort to a perturbational technique. For 
small values of the parameter e = yJrnjM = rQ^JmJl, 
where tq = \/ 1 /M is the radius of gyration of the device, 
the change in angular velocity upon collision with a gas 
particle is small and we can apply a Kramers-Moyal type 
of expansion. Note that this expansion is not uniform 
when performed in terms of the probability distribution, 
but has been found to converge very well when performed 
at the level of the moments. Here we only present the 
final results of this procedure, for more details see [S]. 

Up to first order in e the average angular speed of the 
motor obeys the following equation: 




+ 0{e') 



(4) 



To lowest order in e we recognize a linear drag law, 
Idt{Lu) = —"f{u}), featuring a frictional torque which is 
proportional to the average rotational speed. The pro- 
portionality factor 7 — 7^ is equal to the sum of the 
friction coefficients 7,: contributed by each of the engine 
parts. From Eq. one finds the following explicit ex- 
pression for these friction coefficients: 



7i = mpiVT, 



dU 



(5) 



Note that at this order of the perturbation, no system- 
atic steady state motion appears (w)"* = 0. This indi- 
cates that the "rectification of the fluctuations" leading 
to systematic motion appears at the level of nonlinear 
and non-Gaussian effects. 

At the next order in e, the equation for the first mo- 
ment is coupled to the second moment, whose evaluation 
is thus needed to close the equation. Restricting our- 
selves to the steady state, one finds, not surprisingly, 
that (to lowest order) the average kinetic energy of the 
motor is given by the usual expression for equiparti- 
tion, = ifcTcff, but at an effective temperature 

Toff. The latter is found to be equal to the weighted 
geometric mean of the temperatures in the reservoirs: 

Tcff = {Y.^l^Ti)/{T.^l^)■ 

Combined with Eq. Q we conclude that (up to first 
order e), the engine will develop an average steady state 
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TABLE I: Properties of the 3-d device with a shape optimized for maximum rotational frequency and coohng power, respectively 
(see Fig.[l]). Each unit of the device is cylindrical with parallel surfaces of area vri?^, R — 3nm, height h = 3nm. We assume 
following values of the parameters: density 1350 kg/m'' (typical for proteins); total mass M = 2.29 x 10~^^ kg; mass ratio 
m/M — 1.3 X 10"* (m being the mass of a water molecule); T = 300 K; temperature gradient for the motor: AT = O.IK. 





(^) (Hz) 


7 (10-^* Nms) / (lO"^'' kgm^) 


Qi^a/r (10'^ J/ (Nms)) 


Tiin, (10-21 Nm) 


Q-r (10-15 J/s) 


Motor 


2180 


0.90 1.26 


6.53 


1.17 


1.92 


Heat pump 


1470 


4.55 2.22 


4.41 


4.01 


4.42 



angular velocity given by 



(6) 



We proceed to discuss this first central result of our 
paper. The angular velocity is obviously zero at equilib- 
rium, Ti = Tgff, and also when § dUrf^ = 0, in agreement 
with the fact that the object then loses its chirality. As 
far as maximizing rotational frequency is concerned, a 
numerical procedure was employed to identify the op- 
timum configuration (shape plus axis implantation) by 
deforming the contours in both compartments. This re- 
sulted in the spiral shape depicted in Fig.jljb), with one 
engine part the enantiomorph of the other one. The 
same shape remains optimal, but appearing as the ba- 
sis of a cylindrical object as depicted in Fig.jlja), when 
turning to the case of dimension 3. The dependence of 
the rotational speed on the shape and axis implantation 
is very intricate. In fact even the direction of the net 
rotation is not at all obvious. In Fig. [2] we reproduce, 
for a specific triangular motor element (and its enan- 
tiomorph) , the lines of equal amplitude for the rotational 
frequency as a function of the implantation of the ro- 
tation axis. The rotation is clockwise/counterclockwise 
in the dark/light shaded areas respectively. In view of 
the technological interest of this result, we include the 
corresponding properties of such engine calculated un- 




FIG. 2: (Color online) Equal-amplitude lines of the rotational 
frequency {uj) as a function of the location of the rotational 
axis, for two triangular shapes (apex angle 30° and 45°). The 
(a;, y) coordinates represent the location of the axis with re- 
spect to the center of mass of the unit. Thick lines correspond 
to {ijj) = 0, and dots with maximum frequency. 



der physically realistic conditions, in Table|T] Rotational 
frequencies in the kHz regime are obtained for a temper- 
ature gradient of .1 K. Such a sustained average rotation 
will dominate over the thermal Brownian motion, with 
typical frequency ^JkT/I = 1.81 x 10^ Hz, on a time scale 
of seconds or more. 

We next turn to the analysis of the heat pump func- 
tion. In the following, we will focus only on the linear 
response property, which can directly be obtained by in- 
voking Onsager symmetry. To do so, one needs to write 
the result Eq. (|6| in the framework of linear irreversible 
thermodynamics [6 . One identifies the fiux Ji = {uj) 
and the thermodynamic force X2 — I/T2 — l/Ti. For 
a small temperature difference AT, Ti = T + AT/2, 
T2 = T — AT/2, a linear relation between flux Ji and 
force X2 = AT/T^ is observed, namely: Ji = L12X2. 
The value of the coefficient L12 is found from Eq. ^ 
(for simplicity considering again enantiomorphs) : 



L12 — 



2fcT2 



P\P2 



Ivt (pi + P2Y § dlrl 



(7) 



Following Onsager symmetry 7 , there is a mirror rela- 
tion J2 = L21X1 with an identical proportionality co- 
efficient L21 = T12, while J2 is the flux associated to 
the temperature gradient X2, i.e., it is a heat flux Qi^2 
(from reservoir 1 to reservoir 2), and Xi is the thermody- 
namic force associated with the rotation, namely a me- 
chanical torque divided by the temperature of the sys- 
tem, Xi = T/T. The relation J2 = T2i-'^i, with Eq. (0), 
implies that the heat flux Qi^2 is given by 



^^ 2fcT p,p2 fdlr^, ^ 
Ivt {pi+ P2Y § dlrl 



(8) 



This is the second basic result of this letter. Note that 
the direction of heat transfer depends on the direction of 
the torque, in such a way that it activates an opposing 
Brownian motor in agreement with Le Chatelier's prin- 
ciple [B]. For example, considering § dlrf^ > 0, the motor 
rotates clockwise (w) < for AT < (Ti < T2). The 
application of a positive torque F > 0, inducing counter- 
clockwise rotation, produces an energy flux Qi^2 > 0, 
tending to activate the clockwise Brownian motor. 
The Onsager coefficients L21 and L12 are the off- 
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FIG. 3: (Color online) The cooling power as a function of 
(a) the radius of the device for heights h — 4,6, 10 nm (half 
limiting torque, higher curve corresponds with lower h), (b) 
the applied torque F for height h = 6nm and given radius 
R, (c) the input power Pin for the same dimensions, and (d) 
the location of the rotation axis with respect to the center of 
mass. The shape of the device is optimal for (a,b,c), while for 
(d) it corresponds with the 45° configuration of Fig. [5] Other 
properties as for Table [l] 

diagonal elements of the 2x2 linear response matrix: 

Ji — LiiXi + L12X2, Lii = T/7, 

J2 ^ L21X1 + L22X2, L22 = lil2kT^/{lI)- (9) 

The diagonal elements Ln, the rotational mobility, and 
L22, the thermal conductivity, can again be calculated 
from the above perturbational method [5], or from gen- 
eral arguments based on Langevin theory [3] . As we pro- 
ceed to show, these terms, associated with Joule heating 
and heat conduction, specify the domain in which the 
heat pump can operate as a cooling device. 

An external torque can induce a cooling flux, by pump- 
ing heat out of one reservoir (into the other). This ef- 
fect is offset by a dissipative contribution in both reser- 
voirs, resulting from frictional heating. The linear re- 
sponse term LuXi expresses that work is performed 
on the pump upon application of an external torque F, 
which leads to a power input r^7i/7^ in each reservoir 
i. The dependency ensures that the cooling flux, 
which is proportional to F, will dominate for torques be- 
low a certain Fum, as depicted in Fig.jsjjb) for a con- 
crete realization of the device. The formal condition 
|i2i-'^i| > r^7i/7^ enables us to quantify the limiting 
torque as Fiim = 7^|Li2|/(7ir), which is remarkably 
scale-independent. Maximum net cooling occurs for half 
the limiting torque Fiini/2. Under this condition, a de- 
vice of a few nanometer thickness is capable of a net 
rate of femtojoules per second, cf. Fig.jsja) and Table|l] 
We also note that under this condition, half of the in- 
put power Pin is used for cooling, yielding a coefficient of 
performance rj = Qnct/^in = 0.5. For lower values of the 
torque, F < Fii,ii/2, a higher performance 77 is feasible 



even though, as explained above, the net cooling flux is 
no longer maximal [see also Fig.jsjc)]. 

We finally turn to the issue of thermal conductivity. 
Suppose that a temperature gradient develops under the 
application of an external torque, cooling one reservoir 
and heating the other. The heat pump, being in contact 
with reservoirs of unequal temperature, will then conduct 
heat against the cooling flow. Eq. ([9| tells us that this 
heat flow, J2 — L22X2, has the form of a Fourier law, 
kAT, with conductivity k = L22/T'^- An upper limit for 
the relative gradient AT/T emerges when the conductive 
flow L22X2 equals the cooling power |L2i-'^i|- In explicit 
terms and at maximum performance, F = Fii,ii/2, the 
temperature gradient is bounded by 

AT^Trm p2 {§dlr^^ 
T 8 / {jdlr^y 

The fact that this term is proportional with m/M in- 
dicates that the device may be better suited to transfer 
heat than to create a direct temperature difference. 

In conclusion, a chiral molecular device, operating as 
a heat engine or heat pump, appears to be technologi- 
cally feasible. It remains to be seen whether the simpli- 
fications assumed in this exact theoretical analysis (ideal 
gas reservoirs, frictionless rotation axis and lowest order 
approximation in m/M) lead to a realistic estimation of 
the performance, and whether alternative constructions 
(non-rigid coupling between units, vibrational instead of 
rotational units) offer an even better perspective. 
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